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Abstract 

The low temperature correlation length and the static structure factor are 
computed for bilayer antiferromagnets such as YBa2Cu306- It is shown that 
energy integrated two-axis scan in neutron scattering experiments can be 
meaningfully interpreted to extract correlation length in such bilayer antifer- 
romagnets despite intensity modulation in the momentum perpendicular to 
the layers. Thus, precise measurements of the spin-spin correlation length 
can be performed in the future, and the theoretical predictions can be tested. 
It is also shown how the same correlation length can be measured in nuclear 

magnetic resonance experiments. 
PACS: 75.10.Jm, 74.72.Bk 
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Nearly thirty five years ago Birgeneau, Skalyo and Shirane |IJ invented a neutron scatter- 
ing method, known as the energy integrated two-axis scan (ETAS), whereby the spin-spin 
correlation lengths of layered magnets can be measured without explicitly measuring the 
dynamic structure factor for all energy and momentum transfers and then integrating over 
the energy, known as the three axis scan. Except for special circumstances, three axis scan 
to measure correlation length is difficult. 

The underlying principles of ETAS are an analog energy integration performed by using 
a special geometry and the reliance on the dominance of low energy scattering — the method 
is especially effective when the critical scattering dominates. In recent years, ETAS has 
served us well in unravelling the magnetic fluctuations of the parent compound of one of the 
high temperature superconductors, namely the La 2 Cu0 4 0. The experimental results are 
in excellent agreement with our theoretical understanding of these two-dimensional magnets 
0. Unfortunately, the method does not appear to be readily extendable to a wide class 
of high temperature superconductors with close magnetic bilayers or triple layers within 
the unit cell; a particularly important example is YF^CusOg, which has a close pair of 
magnetic planes within the unit cell, but it is otherwise a square lattice spin S = 1/2 
Heisenberg antiferromagnet. The reason for the difficulty of using ETAS is an intensity 
modulation @] with the momentum transfer perpendicular to the planes. 

It is the purpose of the present paper to (a) derive the low temperature properties of 
such bilayer antiferromagnets and (b) to show how ETAS can be extended to such bilayer 
antiferromagnets. It is hoped that future experiments using ETAS will be able to provide us 
with precise measurements of the antiferromagnetic correlation lengths in these materials, 
which should be of great help in understanding the magnetic flutuations of these bilayer su- 
perconductors. We also briefly remark how the same correlation length can be also measured 
by measuring the Cu relaxation rate in nuclear magnetic resonance (NMR) experiments. 

The Heisenberg model for a spin-5 1 , nearest neighbor, square lattice, bilayer antiferro- 
magnet is 



2 



H= J,, £ Sj p) -Sf + J L £sf ) .S? ) . (1) 

(*j)»p * 

The sum in the first term is over the nearest neighbor pairs in each plane, where the plane 
index p takes two values 1 and 2. The second term represents the coupling between the 
planes. The exchange constants Jy and J± are both positive. 

The low energy, long wavelength properties of the two dimensional Heisenberg model 
is well described by the quantum 0(3) nonlinear a-model 0. Here, we shall consider its 
generalization to coupled bilayers. The Euclidean action for this system can be written 
down on general symmetry grounds ||, but it can also be derived from a (1/5*) expansion 
Hi. The action is 
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where o-W is the component of the staggered field in the direction of the staggered 
magnetic field h. Here, all the coupling constants and dimensional variables have been 
scaled to their dimensionless forms. The label p = 1,2 denotes the spatial directions, and 
p = corresponds to the imaginary time direction. The extent in the imaginary time 
direction is given by u, where Co is the spin-wave velocity. For mathematical purposes, it is 
possible to extend the action to dimensions, d, other than 2 and to the general symmetry 
group 0(N). It is only for d = 2 and N = 3, however, that the model represents a bilayer 
antiferromagnet. Since there are no difficulties in working directly with d = 2 and N = 3, we 
shall use these parameters throughout the paper. For later convenience, we define h = hogo, 
7_l - ^d-, and o t — 

The relations between the u-model parameters and the Heisenberg model parameters are 
known in the large-S" limit B . After correcting minor errors, the correct relations, in the limit 



S - oo, are g = & = hco ^J^ , = and u = fihcoJT^A, 

where 5® = SJ ^+ J± ] the spin-stiffness constants are p° = J\\S 2 and p° = J±S 2 , and the 
spin wave velocity is Cq = ^w2Jj|(4J|| + J±). The quantity a is the lattice constant of the 



Heisenberg model, and the momentum cutoff is given by A = ^JJL to conserve the number 
of degrees of freedom. 

Our calculations show explicitly that the angular momentum coupling between the 
layers is irrelevant for weakly coupled bilayer systems such as YBa2Cu306, although it breaks 
the "Lorentz invariace" and renormalizes the spin-wave velocity. The angular momentum 
coupling term also reduces to a higher gradient term when the number of layers tends to 
infinity. We shall omit this term from our further discussion. 

This bilayer model has a massless mode and a massive mode. We can see this by simply 
expanding the action to quadratic order: 
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where 7I-W is the transverse component of ft^\ In the absence of staggered field ho = 0, 
the symmetric combination is gapless, and the antisymmetric combination has a gap \ 
which will be called the bare dimensionless bilayer gap; the bare dimensional bilayer gap 
is given by hcoK^J~^_. At finite temperatures, the massless mode becomes massive due to 
interactions, which is a topic of this paper. 

We carry out one-loop momentum-shell calculations similar to those of Ref. f| and obtain 
the following renormalization group equations 
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where e is the length rescaling factor. The variable t is the dimensionless temperature 



variable. The dimensionless thickness of the slab in the imaginary time direction u = g/t 
satisfies a simple scaling relation given by (g/t) = (go/to)e~ l . 

The zero temperature flows of the renormalization group equations are shown in Fig.|l|. 
There are two phases, separated by a separatrix between the two unstable fixed points 
(g = Air, g± = oo) and (g = 8tt, g± = 0); the former is the fixed point of the single layer case 
f|. There are two stable fixed points. The ordered-phase fixed point is located at (0,0), 
where both the in- and the inter-plane couplings are infinitely strong. The disordered-phase 
fixed point is located at (00,00), where the system becomes totally disordered. Although 
quantum nonlinear a-model is a very accurate description of the low energy physics in or 
near the ordered phase 0, far into the disordered phase, such a continuum theory can not 
be expected to be valid. 

In the bilayer, 5 = 1/2 Heisenberg model, a quantum phase transition takes place when 



J± > 2.55J|i |nj. In the o"-model, the phase transtion takes place at a critical value of g 
for any g±. At present, we do not have an accurate mapping between the Heisenberg model 
parameters and the a-model parameters allowing us to relate the quantum disordered phases 
of these two models. In the quantum disordered phase, the large-S 1 analysis is not appropriate 
because it relies on the existence of goldstone modes that are not present in this phase. 

The parameters applicable to YBa2Cus06+ x he in the ordered phase, where 7° ~ y- <C 1, 
and g is below the critical value required for the phase transition to the quantum disordered 
phase at T = 0. Therefore, the system is in the renormalized classical regime ||. Also, from 
experiments, it is known that the ground state is an ordered Neel state and that J± ~ 0.1 Jy 

To proceed further, we need an analytical solution of the renormalization group equa- 
tions. We have obtained a good appoximation to the solution based on the following obser- 
vations. In the renormalized classical regime, the bilayer gap Vt± is initially much smaller 
than unity, but increases as V7± oc e 01 "' 1 / 2 , where a (I) < 2, but tends to 2 for large I. We 
can therefore consider two regions, 7_i_ -C 1 in region (I), and ~ 27_|_ in the region (II). 
We solve the renormalization group equations separately in regions (I) and (II), and then 



join the two solutions together to get the final answer. In (I), we replace \/l + 7_i_ by unity 
and obtain 
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where (/i,ii,7jj are the variables at the endpoint of (I). In (II), we use = 2^y± and get 
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The initial condition 7" 1 guarantees that the crossover between (I) and (II) is smooth, 
namely, 7} C 1 and ^r\h — Since we can choose l\ such that 7° <C e" 2 ' 1 <C 1, this 
garuantees that 7^ 1 from Eq. (||). The initial temperature is low enough that both 2; 
and 2^ e_il are m uch larger than unity Following Eq. (|8]), we get ~ ^(1 — |^). At this 
point, the second term in Eq. (0) becomes negligible because -SL coth ^tV /1 +7±) ^ goh_ c -h _ 
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-^z^e' h < 2; therefore ~ 27^_ is also true. 

The two solutions can now be combined together to get 
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The single layer result 
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is trivially recovered when the bilayer gap is much smaller than the inverse correlation length, 
that is, J^l <C e~ l . 



In the limit of large I, the solution has the asymptotic form of 2 = t 2 - — 4-, where 
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This is identical to the low temperature limit of the solution of a classical single-layer a- 
model: = — Thus, we can map the quantum bilayer model to the classical single 



layer model by identifying ^_ — and | = The remaining arguments as to how to 
combine the classical two-loop result, along with the conversion of the pref actor from the 
continuum to the lattice, are the same as those given in Ref. ||. We get 
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with the exact prefactor determined in Ref. |Sj. 

An important result is that the argument of the exponential is twice of that of the single 
layer case. This factor of 2 comes from the mapping between the bilayer model and the 
single layer model. To get a feeling for this factor, it is instructive to consider a classical 
bilayer model given by the action S = £ J d 2 x £? =1 (\d^fl®\ 2 - ha®) + \ - iV 2 ) 
To map it on to the single layer model, we need to define a new field f2, where the transverse 
components are given by n = hi 7 * + 7r ^ 2 ' ) )) an d the longitudinal component is given by 



a = J 1 — | -rr | 2 . This can be accomplished by integrating out the massive fields — n^). 

The result is the new action 5" = / d 2 x (jO^tll 2 — ha) , consistent with our previous result. 

When temperature is low enough such that > 1, the correlation length can be 

written as 

e he ±pM 
« = Sl^Wf , (15) 

where p s (0) is the renormalized spin stiffness constant, which within one-loop approximation 
is given by p s (0) = p° s 1 — ^(1 — \ J^Jl — ■ As discussed in Ref. ||, the formula 
can be used with p s (0) the exact spin stiffness constant of the bilayer Heisenberg model. 

The numerical value of the spin stiffness constant can be obtained from spin-wave theory, 
and it is given by p s (0) = J {l S 2 Z 2 C (S)Z X (S). For S = (1/2), J± = 0.08J|| @, and up to 
order of (1/S), we get Z C (S) = 1.153, Z X (S) = 0.531, and p fl (0) = 0.176 J||. Therefore the 
correlation length is given by 

4^Ps(0) 

^ ~ 0.25ae fc s T , (16) 

where a is the lattice spacing. The exponential temperature dependence is identical to the 
single layer case, but the factor of 2 in the exponent should be noted. 



The static structure factor can be obtained in the same way as in Ref. ||. For the 
bilayer case, it contains two pieces, one for the symmetric combination of the unit vectors 
from each layer and the other for the antisymmetric combination. The notation in the a- 
model is exactly the opposite to the notation in the Heisenberg model, since the unit vector 
field Cl^ in the a-model is the continuum limit of the direction vector of the staggered spin 
operator (— 1)* +P S^ in the Heisenberg model. Our renormalization group equations show 
that in the long wavelength limit both pieces satisfy the following equation as in the single 
layer case || 

21 ( l 0\ c/ J, 



S(k,t )=e M ^j S(e l k,t), (17) 

where t = t(l) is the running coupling constant and k is defined with respect to the antifer- 
romagnetic wavevector (jr/a,7r/a). 

For the symmetric piece, we evaluate the right hand side at I = I* such that we are 
far into the disordered regime, either because the correlation length is of order unity or 
because the rescaled wave vector ke l is sufficiently large. These requirements are satisfied 
by choosing a I* such that £,~ 2 {to)e 21 + k 2 e 21 = 1. When this condition is satisfied, we can 
approximate S s [e l *k, t(l*)] by the Ornstein-Zernicke form 

^'•wi- ^iWr (18) 

Then, 

S.(Mo) = |^/(*), (19) 

where x = fc£ and f(x) = 1+2 1 1 °^ 1 2 +x - . 

The antisymmetric piece can be obtained in a similar manner. When I ~ I*, the bilayer 
gap v^Tl is increased by a factor e l * . At this point, it is safe to assume the Ornstein-Zernicke 
form 

^MC))* g.%) . (20) 



where £j_ is the length scale associated with the bilayer gap. Approximately, we have £j_ ~ 
e ~ P V7±/A ~ \J(1 — f^)7i/A. In terms of the initial coupling constants, we get 

S a (Mo)^|^p. (21) 
The symmetric piece is dominant in the long- wavelength limit because 
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In ETAS, the wavevector of the incoming neutron is fixed, while the outgoing neutrons 
in a direction perpendicular to the layers are collected, regardless of their energies. The 
transferred wavevector is given by q = — q^. Its in-plane component qy is a constant, 
Q|| = — while its perpendicular component is a variable, q± = qj — qi±. For qy near 
reciprocal lattice vectors, the form factor is approximately a constant, and the intensity /(qy) 
is proportional to / °° dqfS 3D (q, ^(qf — if)), where S' 3D (q, uj) is the 3D dynamic structure 
factor. It is related to the 2D dynamic structure factors by <S 3D (q, a;) = sin 2 ( 2 ^)S'f D (qy, u>) + 
cos 2 ( 2 ^ : )S' s 2D (q||, uj), where h is the distance between the two layers. The quantity Sl B (q\\,uj) 
corresponds to the antisymmetric spin combination with respect to the layers, symmetric 
in the a-model sense, and ^^(qy, uj) to the corresponding symmetric spin combination, 
antisymmetric in the cr-model sense. In experiments one probes the region qy ps G, where 
G is the nearest antiferromagnetic reciprocal lattice vector. Defining k = qy — G, we can 
rewrite the intensity in terms of the cr-mode 

1 structure factors, recalling that the definitions of the symmetric and the antisymmetric 
combinations get switched in going from the spin picture to the a-model picture. The 
intensity contains two pieces, I(k) = / s (k) + I a (k), where 

J s (k) ~ f dq f sin 2 (g/ - 2 qi±)h Sr°(K ^(«J - &), (23) 
J a (k) ~ dq f cos 2 (g/ 7 ?±) ^ r (k, ^{q) - g 2 )). (24) 

The gj_-modulation is unimportant in the critical region. The reason is that Sg Da (k, uj) 
is dominated by the critical fluctuations near uj = 0, where both sin 2 ^ 9f an( j 
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= ^Qf are essentially constants. Therefore, we can pull these factors out of the in- 
tegal and obtain the intensity approximately proportional to the static structure factor, 
J s (k) ~ ^- sin 2 ( fe 7' x) /i) f™ E . duoS^ik, u) ~ S B (k,t ), where we have reverted to the 
previous notation by dropping the superscripts. The quantity is the incident neutron en- 
ergy. For the antisymetric piece, we get an upper bound by neglecting the factor cos 2 ( 2 ^ 1 ), 
which is J oo c%S' 2D<T (k, ^(qj - qf)) ~ ^^(k, t ), where (q f ) is some average of the 
wavevector. Because S s (0,t ) ^> S^O,^), the intensity is dominated by the contribution 
from the symmetric piece. Therefore ETAS for a bilayer should yield information about the 
symmetric piece, hence the correlation length. The contribution of the antisymmetric piece 
should result in a small broad background. 

The NMR relaxation rate for the in-plane Cu site will be given by |Tl| ^- oc T§ £, because 
the symmetric structure factor will dominate when the correlation length £ is large. Thus, 
the formula for the correlation length, Eq. (|T5|), can also be tested in NMR experiments, as 
for the single layer La2Cu04 [fL2| . 

In conclusion, we have computed the low temperature properties of bilayer antiferromag- 
nets using a quantum nonlinear a-model approach. Explicit results were given for the static 
structure factor and the correlation length. The correlation length diverges exponentially as 
T — > 0, but the striking feature is that the argument of the exponential function is a factor 
of 2 larger than the corresponding single layer case. We have also shown that ETAS and 
NMR can be used to test our theoretical predictions concerning the correlation length. 

We thank O. Syljuasen for discussions. This work was supported by the National Science 
Foundation, Grant. No. DMR-9531575. 
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